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Abstract 

We have applied a double decoupled localized level Anderson-Newns Hamiltonian 
to the analysis of surface effects upon the ionized fraction 1Z+ of sputtered atoms 
from a metal surface. Electronic excitations, induced in the conduction band by 
the transient formation of quasi molecular systems, between substrate and emitted 
atoms, in the collision cascade generated by the primary incident beam, have been 
explicitly included into an instantaneous transition matrix peaked at the Fermi level 
of the material. The interaction dynamics seem to take place over two different 
time scales, one related to sputtered atom trajectories and the other to recoiled 
substrate particles. Finite temperature calculations have suggested, at very low 
ejection energies, a power law dependence of the final charge state of the sputtered 
beam on its detected velocity. This result is in agreement, in the zero temperature 
limit, with some previously published papers and its validity has been compared to 
other theoretical outcomes and tested on SIMS data. 

pacs: 79.20.-m, 71.23.An 

keywords: ion-solid interactions; ion emission; metallic surfaces; models of non 
equilibrium phenomena; scattering and sputtering. 



1 Introduction 

Resonant electron transfer during sputtering of atoms and molecules from 
metal surfaces has been a very interesting subject in relation to secondary 
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ion mass spectroscopy (SIMS) where the knowledge of ionization and neu- 
tralization rates is essential for many quantitative analyses. The theory of 
such processes should encompass a detailed description of charge exchanges 
in the collision cascade to the outgoing atom or atomic cluster. In fact, to 
make this complicated problem tractable, a basic model has been formulated 
in terms of a single orbital parametrically time-dependent Anderson-Newns 
Hamiltonian [1]-[16], in which the effects of the bombardment, including the 
substrate motion, the electronic excitations generated and the physical and 
chemical consequences of the void produced by the ejected particle have been 
neglected. Many experimental results have been reproduced in this way [17]- 
[19] which is common to both SIMS and atomic surface scattering. Yet the 
need for a more comprehensive theory of surface influences in sputtering has 
led to a generalized approach [16,20-22] in which the substrate dynamics have 
been modeled with locally time-dependent perturbations created during the 
collision and acting within the first layers of the surface region. In a first simpli- 
fied treatment [20] a one-body time dependent scattering potential, localized 
in the vicinity of the emission site, was added to the basic Hamiltonian and 
direct (i.e. non mediated by the atomic localized level) electronic excitations, 
created in the conduction band by this new term, were explicitly considered in 
the evaluation of the ionized fraction. Although the average lifetime of these 
excitations is believed to be rather short on the time-scale of the emission 
trajectory [23], their contribution to the final charge state of secondary ions 
becomes important at very low emission energies where the ionization prob- 
ability calculated from the basic theory exhibits an exponential dependence 
on the inverse of the average outgoing velocity of sputtered particles [6,10]. 
In a recent paper [22], an exact formal solution of the two potential problem 
formulated in Ref.[20] has been obtained confirming more rigorously that, in 
the velocity regime typically observed in many SIMS experiments, the basic 
hopping mechanism may be negligible with respect to surface induced excita- 
tions, at least if they can be represented in the form of a one-body scattering 
potential. In the present article, we have proposed a more accurate form of the 
sputtering potential combining the above mentioned impurity scattering con- 
tribution with a hopping resonant interaction due to a second discrete state 
embedded in the conduction band (Sec. 2). Both perturbations are essentially 
generated in the quasi-molecules (QM) formed transiently in the collision cas- 
cade, between substrate and emitted atoms, and the latter has been already 
used in a different context concerning the analysis of electronic emission from a 
bombarded metal substrate [16,21]. The reference Hamiltonian is quadratic in 
electron creation and annihilation operators, which would make the problem 
trivially soluble in the absence of time-dependent interactions. Nevertheless, in 
the present treatment, only a formal expression can be given to the exact evo- 
lution of the localized orbital on the ejected particle. We have performed this 
step by using retarded Green's functions thus parameterizing the evolution 
of the atomic state in terms of Feynmann diagrams (Sec. 3). Assumptions of 
the slow, non adiabatic dynamics, i.e. Brako and Newns' approach [10,11,24], 
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plus an approximation due to the transient nature of the QM, has allowed 
us to obtain a manageable expression for multiple interaction transition am- 
plitudes involved in the problem (Sec. 4). We have generalized former results 
[20,22], with particular attention to the dependence of the ionization prob- 
ability on the final emission velocity of the outgoing ion, including thermal 
interactions. We have also proposed a comparison of our predictions with pre- 
viously published theoretical calculations [4,6,10,12,20,25] and applied them 
to a well known sputtering experiment [12,26]. 



2 Formulation of the Problem and Generalized Hamiltonian. 

The generalized model is sketched in Fig.l. The metal is treated as a Fermi 
gas specified by a set of continuous states {\k)} and spectrum {ek} with work 
function 0. The Fermi energy has been set to £/ = — in an appropriate scale 
with the vacuum fixed to zero. The multiple index k runs over single particle 
energies and all the degeneracies of each level. In addition, because of impuri- 
ties, adatoms and surface defects, it may include discrete quantum numbers. 
The outward atom is assumed to move along a classical trajectory and have 
a non degenerate valence level e a (t) lying within the metal conduction band. 
The corresponding valence state is indicated by \a(t)). e a may be an affinity 
level (A) as well as an ionization level (I) which give rise, respectively, to 
negative and positive ionization. An electron can tunnel between the atomic 
level and a continuous state \k) with a coupling strength given by the hop- 
ping integrals [1,2] Vk a (t) and V a k(t) = V£ a (t). Particle exchange with the 
continuum dominates in a range of the order of the virtual broadening of the 
localized state 

A a (t) = 7i^2\V ka (t)\ 2 5[e k - e a (t)] , (1) 
k 

which is proportional to the first order transition rate of all occupation chan- 
nels \k) — > | a), i.e Fermi's golden rule [3]. Anderson- Newns' potential causes 
an indirect production of electron-hole pairs by processes [11,16,20,21,27] of 
the kind shown in Fig. 2a. Other excitation channels are induced by tran- 
sient perturbations generated during the collision process. In relation with 
the schematic of Fig.3, we observe that the majority of low energy atoms are 
emitted because of indirect binary collisions with substrate atoms that have 
received a sudden impulse in the collision cascade and the number of ejections 
caused by direct collision with primary ions may be neglected. For this rea- 
son, we consider secondary emitted atoms to form quasi molecular systems 
with recoiling substrate atoms moving within the surface region. With this 
mechanism an electron may be trapped, as illustrated in Fig. 4, into a transient 
molecular orbital (MO) which may be promoted to the conduction band as a 
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result of the increasing nuclear charge in the QM with decreasing internuclear 
distance and the increase in kinetic energy of the electron with its increased 
localization. Thus the hopping interaction of this state with the continuum 
can be very important at low incident energies when the QM has a longer av- 
erage lifetime. In addition, since translational invariance is broken down both 
by the substrate and secondary atomic motion, many scattering transitions 
occur in the conduction band together with electronic emissions. 

In this paper we have limited our investigation to the influence of a single 
promoted state \m(t)), with energy e m (t) in resonance with the conduction 
band, interacting through a second hopping potential, specified by matrix 
elements Vk m (t) and V m k(t) = V£ m (t), with the continuum. This MO may 
be either a metastable configuration of the outgoing atomic state or distinct 
from the atomic state itself, depending on the energy position of the localized 
levels of the two atoms involved in its formation. 

Let us denote by X (t) and Y (t) the classical position vectors of the substrate 
and emitted atoms in the QM, respectively, as measured from a reference frame 
placed at the surface. Because of the slow motion of the atomic cores on the 

— # — # 

electronic time-scale, X and Y introduce two large characteristic times, tx and 
Ty, over which many electronic transitions occur. Moreover, since the substrate 
recoling atom in the QM, see Fig. 3, is more energetic than the emitted atom 
and follows a complicated path inside the metal, we may take Ty >> tx- 

In order to construct the second quantization Hamiltonian of the problem 
we consider a single electron of the system, with kinetic energy ex = p 2 /2, 
p being its momentum operator, in the surface barrier of the metal v s (z), 
with z its orthogonal coordinate operator, and in the screened Coulomb field 
v c T — Y(t) +v c r — X(t) , of the atomic cores in the QM, with Fits position 
operator. Such a particle is also subject, at large distance form the surface, to 
the potential of its image charge and of those of the screened nuclei at (y, X^j . 

Denoting with vi r; Y (t) , X (t) the sum of these three contributions, we can 
write the first quantization Hamiltonian of the particle as 



h (t) = e K + v s (z) + v c r-Y(t) 



+ V c 



X(t) 



+ Vi 



r;Y(t),X(t) 



(2) 



where atomic units have been used. Next, we define the continuous con- 
duction states {\k)} as eigenstates of the unperturbed surface Hamiltonian 
h s = ex + v s (z), with spectrum Because of the surface barrier, the cor- 
responding wave functions (r\k) drop to zero outside the surface, i.e. at z > 0, 
after some atomic distances. The adiabatic state \a(t)) is obtained by trans- 
lating the stationary state |a°), for an electron localized on an atomic nucleus 
fixed at the surface, to the position of the emitted atom. Hence it diagonalizes 
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instantaneously the Hamiltonian h a (t) = ex + v c r — Y (t) , with a time inde- 
pendent eigenvalue e° a = —I (or e° a — A). The atomic wave function depends 
on f — Y, i.e. (f]a (t)} — (r — Y (t) |a ^. As for \m (t)), it can be calculated by 
solving the eigenvalue equation for the QM Hamiltonian 



h m (t) = e K + v c i — Y(t) 



+ v c 



X(t) 



viz. h m (t) \m (t)) = e° m (t) \m (t)), through the LCAO method. 

The set of normalized states {\k) , \a (t)) , \m (t))} is well defined and can be 
used as a truncated basis to describe the problem, once the overlap contribu- 
tions (a (t) \m (t)), (a (t) \k) and (m (t) \k) have been neglected for all t. Thus, 
we can introduce the electron destruction (cj)j e ^ amk y and creation operators 

, in Schrodinger's picture, satisfying algebraic relations 



j&{a,m,k} 



(3) 



with |...| the anticommutation operator and rij = CjCj electron number op- 
erators. These may be chosen to be approximately time-independent by an 
appropriate phase transformation on \a(t)) and \m(t)) and the many-body 
system Hamiltonian [28] is 



H (t) = e a it) n a + e m (t) n m + J2 e kn k 

k 

+ T,{ Vka (*) C l° a + h ' C -} + E{^™ (*) C l°m + h-C. J 
+ J2 V kk' (*) c{c k >. 



(4) 



k,k' 



Here the atomic and molecular unperturbed energies, Eq and e™, are shifted 
in presence of the interactions to 



e a (t) =e° a + (a (t) \ {u c [r - X (t)} + v s (z) + v z [r; Y (t) ,X (t)} } \a (t)) 



and 



e m (t)~e° m (t) + {m (t)\v s (z) \m(t)), 



(5) 



(6) 



respectively. The hopping terms take the form 
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V kj (t)~(k\{v c [?-X(t)]+v c [?-Y(t)]}\j(t))^ with je{a, 



and the scattering potential can be written as 



V kk ,(t)~{k\\v c r-X(t) 



+ v c 



f-Y(t) \\k'). 



(8) 



In Eqs.(6-8) we have omitted image potential contributions which act sensibly 
when the emitted atom is at large distance from the surface where the contin- 
uous wave functions vanish and the MO does not contribute. Non Anderson- 
Newns terms enhance the number of electron-hole excitations by processes 
similar to those reported in Figs. 2b and 2c. 

Moreover in the reference Hamiltonian we have not allowed direct transitions 
between the two localized states, regarding terms like 



V am (t) = (a(t)\v s (z)\m(t)) = V: m (t) 



as small perturbations in the energy scale of the hopping interactions. We 
point out that the introduction of the QM contribute in the determination 
of the quantities e a and V ka that, unlike the bare model [3,29-31], depend 
parametrically on both the outgoing and recoiled atomic trajectories. 

Brako and Newns' slowness approximation [10,11,24] can be formulated ex- 
tending its application to the problem of an electron gas probed by two local- 
ized sources. We assume that changes of matrix elements ({j\h j>e{k am}' 
over the characteristic times of the atomic motion in the QM, are negligible. 
More rigorously we have quasi- stationary dynamics such that 



(j\h\f) 



(J\h\f) 



[Y(t),X(t)\ 



[Y(t),X(t)} 



^(j\h\f) 



~(i\h\f) 



[Y{t+TY),X(t)] 



[Y(t),X{t+T X )\ 



(9) 



As we have already observed surface effects rapidly dissipate with respect to 
Anderson-Newns hopping processes, and a localized molecular electron has a 
short average life-time in the Ty-scale. Then, for times larger than ry, all the 
terms corresponding to inner excitations may be neglected. However while the 
QM is present many electronic exchanges, see Figs. 2b and 2c, occur that may 
influence the final charge state of the sputtered atom. 



6 



A very important topic is concerned with the question weather the two discrete 
states are so much linked that \a) may transform into \m), in the first step of 
the emission process, during the action of the short range atom-atom forces, 
or if they can be considered distinct and non overlapping. In our opinion, the 
former case can occur when the energies e a and e m are so close, at the time 

— * — * 

of ejection, viz. t* = 0, that one of the atomic states, localized at X and Y 
respectively, in which |m(0)) may be decomposed, in the LCAO scheme, is 
| a (0)) itself. Anyway, even in this eventuality, we can continue to use the model 
Hamiltonian of Eq.(4), assigning to the state \m) the role of a discrete atomic 
state localized on the substrate atom of the QM and transiently forming a MO 
with the state \a) in the time interval < t < tx- Apart from this conceptual 
correction no significant qualitative or quantitative changes affect the ionized 
fraction in the emitted beam, at least within the limits of our approximations. 
Hence, in the following, we shall consider \a) to be independent from \m), 
regarding a more accurate study of correlations between the two discrete states 
to future works. 



3 Exact Diagrammatic Solution 

The purpose of our analysis is to find out the fraction of atoms detected at 
the end of the process with a given charge state. In our model this calculation 
is straightforward once we have determined c a in the final system configura- 
tion, when the particle is quite far from the sample and the perturbations 
are switched off, i.e. c a (oo). The Hamiltonian (4) describes atomic emission 
only after the collision has taken place hence the initial condition should be 
referred to the time when the atom is at the surface. For analytical simplicity 
we prefer to place the initial time in the remote past and switch on adiabat- 
ically the perturbations as in an ion-surface scattering problem. It has been 
argued that this analytical continuation on negative times does not infer the 
treatment for many systems, such as Na/W(100) [11]. We therefore assume 
that interaction terms obey the following 



Vkk' 


( 


— oo) 


— Vkk' 


(oo) 


= 


v ka 


( 


— oo) 


= V ka 


;<») 


= 


Vkm 


( 


— oo) 


= V km 


(oo) 


= 



The neutral fraction of emitted particles at a temperature T is defined as 
the thermal average of the atomic electron number operator, evolved to the 
remote future according to Heisenberg's scheme, in the initial equilibrium 
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configuration, or 



Vq (T) = (n a (oo)) . (10) 

Here, spin degrees of freedom have been ignored as Coulomb intra-atomic 
correlation terms [24,32-35] do not explicitly appear in the reference Hamilto- 
nian. This approximation limits the applicability of the model to very simple 
ejected beams such as alkalis [10], and, even in this case, there are non triv- 
ial problems at least for the atomic orbital. In fact, while the emitted atom 
moves off surface, the state \a) increases its localization, the hopping mecha- 
nism becoming weaker, and approaches the atomic limit, when the exchange 
of electron with the surface becomes a spin-flip process without charge fluctu- 
ation [24], which cannot be described by a spinless theory. On the other hand, 
the state \m) acts sensibly on the system only when the emitted atom is in 
the vicinity of the surface so that its correlation energy may be taken to be 
infinite, in the scale of hopping interaction strengths, avoiding double occupa- 
tion of the level. Some significative efforts in a better comprehension of spin 
correlations have been made, in the last decade, on the bare Anderson-Newns 
model, see for instance Refs. [33,35]. However further studies are needed to 
discuss the influence of the two-body spin-spin potential in the framework of 
the present approach. Thus, in this analytical derivation, the spin degeneracy 
can be taken into account only by addiction of a 2, premultiplying Eq.(10), at 
the end of calculations. Complementarily to Eq.(10) the ionized fraction can 
be put in the form 

TZ + (T) = 1 — V (T) = (c a (oo) 4 (oo)) . (11) 

Our evaluation of (10) and (11) are based on the equations of motion method 
[7]. Electron operators obey the following set of coupled equations 

i±c a it) = e a it) c a (t) + E V ak (t) c k (t) 

k 

ij;C m (t) = S m (t) C m (t) + E Vmk (t) C k (t) 

k , (12) 

if t c k (t) = e k c k (t) + E Vkk' (t) c k ' (t) 
k' 

+V ka (t) C a (t) + V km (t) C m (t) 

that can be easily written down by using the algebraic relations of Eq.(3) 
extended, in Heisenberg's picture, at all equal times. It can be noted that 
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the discrete states \a), \m) are coupled only through indirect processes, as in 
Fig.2d. 

Eqs.(12a-c) can be combined and iterated to give a formal result in terms 
of time dependent Feynmann graphs: the integral form of (12b) is obtained 
by introducing the unperturbed retarded Green's function, or empty band 
propagator, of an electron in a state \j), i.e. 

i 

—i j drSj (t) 

G °+(t,t') = -t9(t-t')e * 



as solution of the singular equation 

j^- £ ,(*)j G o+ (t,t') = 5(t-t>). 

Then, we have 

oo 

c m (t) = iG°+ (t, -oo) c m + J2 J dt'G°+ (t, t>) V mk it') c k (f) . 

k — oo 

Substituting into Eq.(12c), the new equation for c k (t) reads 

oo oo 

c k ( t )=tG° k + (t,-oo)c k + Y: J dtf J dt»G° k + (t,t')W™, (t',t")c k > (t") 

^ — oo — oo 
oo oo 

+i J dt'G k + (t,t')V km (t / )G° rr t(t',-oo)c m + J dt'G° k + (t,t')V ka (t')c a (t'), 



(13) 



where surface-induced excitations of Figs. 2b and 2c have been enclosed into 
the non instantaneous potential 

(t, *') = V kk , (t) 5 (t - t>) + V km (t) G°+ (t, t') V mk , (*') . (14) 
Eq.(13) may be substituted into Eq.(12a) to obtain 

{i-- £ „(0}c a (0- / dt'i; a (t,t')c a (t>) = iY / v ak (t)G k + (t,-oc)c k 
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oo oo 



+ £ /df / dt»V ak (t)G° k + (t,t')W^(t',t»)c k ,(t») 

oo 

+*E / dt'V ak (t)G° k + (t,t')V km (t')G°+(t',-oo)c m , (15) 

k 

" — OO 

in which S° (t,t') is the retarded self-energy of the state \a) in the absence of 
the QM, 

K (t, = E ^ (*) (* - v ka (0 . (i6) 



This term is peculiar to the Anderson-Newns Hamiltonian and responsible for 
broadening and shift of the atomic level. Surface induced excitations act to 
renormalize S° and G® + ; after an infinite number of iterations of Eq.(13) into 
Eq.(15), we can specify an inner retarded propagator, 

(t,t) = G>» (t,t')5 k , k , 

oo oo 

+ E f dt ± [ dt 2 (ff(t,t 1 )W& fat*) G]ft faf), 

— oo — oo 

representing the probability amplitude of a process \k') — > from time t' to 
time t, with multiple empty band scattering from WJJj:,, and we can rewrite 
Eq.(15) as 



[ij t -e a (t)}c a (t)- J dt'Z a (t : t')c a (t')=iJ2V ak (t)G™t(t,-oo)c k , 

— oo kk' 

oo 

+tJ2 J dt'V ak (t)G^(t, t')V k , m {t')G°+{t\-oo)c m (17) 

— oo 

The dressed self-energy graph is then given the form 



S a (t, f ) = S° (t, if) + TP a (t, ?) = E V ak (t) G™t (f , W« (0 , 

kk' 

(18) 

where is the renormalization term 



oo oo 



(19) 
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and 



(t,t') = W™, (t,f) 



oo oo 




the retarded transition matrix in the QM potential. Integration of (17) gives 
the formal result for the full atomic electron destruction operator 

c a (t) = iG+ (t, -oo) c a + iG+ m (t, -oo) c m + i ^ G+ fc (t, -oo) c k , 

k 



in which the full retarded propagators appear 



oo oo 



G + a it, f) = G°+ {t, t>) + / dt! / (ft 2 G2 + (*,ti)E a (t 1 ,t 2 )G+(t 2 ,* / ) 

— oo — oo 

G + ak (t, if) = E 7 dhGi (t, h) V ak , (tO G™+ fr, 

fc' — oo 



G+ m (f , f ) = E / dtiGi (f , *i) V km (h) G°+ (t u t>) 

k — oo 



Then, the ionization probability is 



TZ+ (T) = G+ (oo, -oo) (ra ) + G+ m (oo, -oo) (n m ) + J2 pafc (oo, -oo) (n fc ) 

(22) 



with (n k ) the unperturbed hole distribution function 

1 



(n k ) 



j(e f -e k ) 



+ 1 



and (n a ), (n m ) the initial vacancy of the atomic and molecular orbitals. Eq.(22) 
generalizes the result of Blandin, Nourtier and Hone [6] on the basic Anderson- 
Newns model obtained on the Keldish contour [36]. 
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4 Approximations on the slow dynamics 



A simplified solution for Eq.(22) can be obtained by repeatedly performing 
the following approximations on each term of the perturbation series defining 
the retarded propagators of Eqs.(21): 

(1) Let us consider the integral 

t 

* -if dr[e k -n(r)] 

X (0 = E J dt 'h (0 e * 

k -oo 

where /& and is a slowly varying function of time, on the electronic scale, 
as specified in Sec. 3, and Q = e a or Q = e^. The oscillating complex 
exponential in the ^'-integral varies sensibly with t' respect to fk so that, 
in the summation over the continuous modes k, it contributes a quasi 
delta function of t — t', at least until the outgoing particle remains in 
the interaction region with the surface. This allows us to hold /& and Q 
constant during integration [10,11] approximating 

with G° + [Q (t)] being the Fourier Transform (FT) of the unperturbed 
retarded Green function G° + (t) in the e = VL (t)-domain, viz. 

G° k + [n (*)] = 7T-- — ^ — - r] -> o+. 

(2) Each term of the kind 

t 

* -if dr[e m {r)-n{r)} 

y m (t)= J dt'f m (t')e «' 

— oo 

where / m contains the inner hopping matrix elements \4 m , V^* m simplifies 
to 

y m (t) = l fm (t) G°+ [Q (t)] . 

In fact for times t < t' < Tx the functions f m ,s m and f2 are slowly varying, 
i-e. / m (£') = / m (i) and 

t 

J dr [e m (r) - SI (r)] [e m (t) - ft (*)] (f - *') , 
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while f m switches off to zero for t' >> tx- We point out that for notational 
convenience we have omitted other quantum indices or time labels on 
which f k and f m may depend. 
(3) The operator 



Q a (t) = J dt'X a (t,t')c a (t'), 



may be treated by replacing c a (f) with iG° a (t 1 , t) c a (t), with G° (t 1 , t) 
the free advanced propagator of the outward atom electron, i.e. 



Q a (t)^i J dt'Ea (t,t')G°- (t',t)c a (t) 



In fact perturbation changes in the full propagator G~ (t', t) are negli- 
gible in the neighborhood of t where E a , as a function of t', reaches its 
maximum. This last is completely equivalent to the semiclassical approx- 
imation (SCA) used in Ref.[35]. 

After some tedious algebraic manipulations, retarded diagrams, introduced 
in Sec. 3, are converted, by means of the approximations introduced above, to 
instantaneous diagrams. In particular the effective QM potential, the T-matrix 
and the renormalized self energy contributions, Eqs.(14),(18) and (20), become 



(*» t') - W Tk> (e, t) = V kV it) + V km (t) (e) V mk (t) 

(t, f ) - T™, (e, t) = (e, t) + £ W£, (e, t) G° k t (e) T™ k , (e, t) 

k" 



k 



(t, t') - Eg (e, t)=E V ak (t) G°+ (e) T™, (e, t) G° k t (e) V k , a (t) 
k,k> 



where T^., (e, t) is a Brako and Newns instantaneous T-matrix [10] and the 
real and imaginary parts of S a (e, t) give respectively the instantaneous shift 
and broadening (see Eq.(l)) of the atomic level in interaction with the QM. 
Eq.(17) is transformed into the integrable form 
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{ i Jt~ £a W -Sa ^}° a (*) 

lim {Y,Vak(t)e 



8 k ,k>-i J dt'e^-^'T™, {e w ,t') 



t 

-i E Vak (t) e-^ J dt'e^-^'T^ [e m (t) , t'} e^^cjj, 



with S a (t) standing for S a [e a (t) , t] and 

(*),*] = E 

(t),*] G°/[5 m (*)]y fc , m (*). 

it' 



(25) 



Eq. (24) can be solved by introducing the semiclassical propagator of an elec- 
tron in the state |a),i.e. 



-i J dr[£ a (r)+S a (r)] 

G% + (t,t') = -iQ(t-t')e f 



Hence the slow solution for the Heisenberg atomic electron operator is 



i 

_ l r -if dr[e a (r)+E a (r)] 

Ca(f)= lim E / " KtM^ 



x 



e i£k ' to c k/ , 



where we have ignored the terms containing information on the initial vacancy, 
referred to the remote past, of the atomic and molecular orbitals. These as- 
sumptions fit the case understudy well. In fact while atoms move away from 
the surface so many exchanges occur with the conduction band that the single 
transition processes \a) — > \a) and \m) \a) may be ignored. From Eq.(ll) 
we have the ionized flux as function of instantaneous terms that can be ap- 
propriately modeled according to the outgoing ion and surface motion 



^ + (T) = E 

k J 



°f » J <M £ fc'- £ a( r )- S a( r )] 

/ ^E e * V ak'(t) 



X 



5 k , k >-i J dt'e^'-^'T^ (e k ,t') 



\n k ) 



(26) 
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Surface effects are included within the instantaneous multiple transition term 
T^ fc acting to renormalize the instantaneous self energy and producing an ad- 
ditional term within the square modulus in the r.h.s. of the previous equation. 
If we switch off the QM potential, letting W™ fe , — > and hence T™ fe , — > 0, we 
reattain the well known result [6,10,13] 

K + (T) = £ 

k 

Our analytical investigations focus on the non adiabatic ionized flux at very 
low emission energies, i.e. 1 -< E -< 20 eV. Non adiabaticity is achieved by con- 
straing the imaginary part of the instantaneous self energy to vanish asymptot- 
ically more rapidly than 1/t [6]. In addition we believe that, within the energy 
range considered, only the electron-hole excitation matrix element in the r.h.s 
of Eq.(26) contributes, in other words continuous electrons are strongly per- 
turbed by W^./ undergoing multiple scattering processes before hopping into 
the atomic level. For mathematical, simplicity we observe that both V km and 
Vkk> are much more sensitive on X, so that we may neglect the F-dependence, 
or assume that QM potentials are effectively localized on the substrate atom. 
The unperturbed substrate may be approximated by a free electron model 
with k representing the eigenvalues of the one-electron wave vector operator, 
i.e. k = k. Indeed k should be a complex vector, taking into account the de- 
caying part of continuous wave functions outside the surface, yet we neglect 
its imaginary part since the QM interactions are peaked in the space region 
where a metal electron can be considered to be in a plane-wave state, there- 
fore this change does not modify substantially our calculations. We can now 
perform a phase transformation on the instantaneous T-matrix of Eq.(23b) 
expressing the ionization probability in terms of a reduced T-matrix localized 
at the surface 

7$ (e k ,, t) = e ^-k')-x { t) T ^ (ejfc/> t) . (28 ) 

Next we assume separability of V ka between k and t dependence [1]-[16], 
Vka(t) = ^jfc| u a(^), where u a is entirely determined by the atomic emission 
trajectory, ignoring the correction due to the substrate atom, and, referring 
to the simplest case, we adopt a rigid trajectory parameterization. We re- 
place Y by its average (y^ = (vt) t so that u a = exp (— |vt), in which v is 
the component of the average outgoing velocity Vx in the orthogonal direc- 
tion to the surface plane [23,14] and (7v) _1 a characteristic time of the order 
of Ty. Moreover, working under resonance conditions of e a and e m with Ef, 
the single particle density of states p as well as the stationary part of the 
Anderson-Newns hopping potential can be fixed at the Fermi surface [6,10] 
and S° becomes a purely imaginary quantity, i.e. S° (t) = —iA a exp (—yvt). 



y i J dT[e k -e a (T)^ a (r)] 

J dte > V ak it) 

-oo 



{n k ) 



(27) 
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Furthermore we consider scattering transitions induced by the reduced inner 
potential of Eq.(23a) to be isotropic, see Eqs.(7), (8), and we use the approx- 
imated relation 



1 



^ e - i£ fc'(*-*')-*fc'-^M7j™ (e k ,t') = 



with kf the Fermi wave vector, obtaining 



sin 


'kf 


X(t) 


kf 


X(t) 





CT £) (e k ,t)6(t-t), 



r dp r sin h 

(T) = Anp A a ( —7- %r / dt— L 



sin 


kf 


X(t)~ 


kf 


X(t) 





re/ fee v ' 



(29) 



in which the correlation term 
SP (t) - —npA a 



sin 




X(t)~ 


*/ 







(30) 



has been omitted. We have thus selected charge exchanges, mediated by the 
Fermi level, of the type reported in Fig.2d. The main part of the process takes 
place over the r^-scale. In the hypothesis ^ 3> t x we can replace the factor 



F(t) = e-i vt exp -e 



a —■yvt 



7v 



within the t-integral in Eq.(29), with its maximum value F max ~ y sa~- m 
addition e a (t) is evaluated at the time t = ln(2A a /7v) /7V, when F (to) = 
F max . Consequently we have the finite temperature generalization of Eq.(16) 
of Ref . [20] (apart from a factor 4 due to a typographical error) 



47T 

U + (T) = — 1 vN T (e a ,e f ), 



(31) 



where 



N T (e 



a,£f) = f 



de 



,p(s f -e) 



+ 1 



00 
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}f 
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kf 


X(t) 





r- rk Ae,t)e-^ 



(32) 
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is the number of holes in thermal equilibrium excited to the Fermi level by 
the effective potential WJIjj./. Our problem has been reduced to a time param- 
eterization of the instantaneous T-matrix. We note that the ionized fraction 
is directly proportional to the atom velocity in contrast to the exponential 
dependence on v _1 obtained at higher kinetic energies and low temperatures 
[6]. This may be taken as an a posteriori justification of the approximation 
that let us neglect the first term in the r.h.s. of Eq.(26). 

A simple analytical formula can be obtained if in W^., we retain only the 
hopping term V km due to the promoted MO. Hence the series (23b) can be 
summed exactly to yield 

_ v kfm (t) G^»v mN (t) 

T *l M) " l-G^(5)E m M) ' (33) 



where S m is the instantaneous retarded self-energy of an electron in \m). We 
observe that X (t) may be replaced with its average value, for < t < Tx, 
when the t-integrand in Eq.(32) is sharply peaked. In the same interval e m 
depends linearly on time, with a slope b depending on its position relative to 
Ef (see Fig.l). Indicating with ti the crossing time, we set 

e m {t) =e f + sgn (e af ) b (t - h) , (34) 



with e a f = s a — Ef. Moreover Vigi can be evaluated at the average X = 



X (t) . Thus we can use the expression 

r r\t i (e, = r-%7 \ — 7T-, (35) 



with A m = Tip Vk f m being the maximum broadening of the promoted MO, 
into Eq.(32) and continue to integrate from — oo to oo for analytical simplicity. 
Employing the FT of (35) in the (e — e a )-domain we find 



— oo 

oo 



+0i 



B(^,/*^i}, ,aa) 



having performed the change of the variable ^ — £ — Ef. Eq.(36) can be inte- 
grated analytically, neglecting contributions of the order of e -/3 l £a/ l, to obtain 
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167T 2 A^ 7 V 2^\e af \ 

56 2 /5sin(2vrA m /6/3) 

(37) 

Such an approximation is referred to the case of a soft promotion, with 
b/2A m ~ 0.1 eV for typical values of \e a f\ of the order of 1 eV and tem- 
peratures up to 1000 K. Moreover when thermal effects become negligible, i.e. 
A m /fo/3 < 1, we obtain [22] the generalization of Eq.(47) of Ref.[20] 

^ ^ , / \ 87i"7vA„, ZsinA;fX\ 2 9 A m i, i 

^ + ^ (e a/ ) - sgn (e a/ ) ^ m j e~ 2 ^M, (38) 

including the case e / > 0. Eqs.(31),(38) are to be compared with the results 
of Brako and Newns (Eq.(34)of Ref.[10]) 



holding in the 1/(3 7t>/7r-limit, and Blandin, Nourtier and Hone (Eq.(61) 
of Ref.[6]) 

11+ - (e af ) - -sgn (e af ) e~^\^ . (40) 



These can be obtained from Eq.(27) under the same assumptions that lead to 
Eq(31). 



5 Comparison of Ionization Probabilities and Conclusions 

In deriving a rough estimation of the fundamental parameters of the problem 
we note that the characteristic energies eo = 'yv/n and e\ = b/2A m are of the 
order of 1/ry and 1/tx respectively. Thus, under the working hypotheses of 
this paper \e a f \ must be larger than both e\ and e while we should also have 
£i ^> Eq. The evaluation of the average outward velocity v requires further 
consideration: it has been well established that when the sputtered ion is in 
the interaction region of the surface its actual velocity is higher than its final 
measured value Vf owing to the attractive potential. In particular this should 
be true in the low emission energy range under study. According to Lang's 
model [15,23], the binding potential of the outgoing atom, of mass m, with 



K + (T) = e(e a/ )-sgn(e o/ ) 



sin kfX 
k f X 
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respect to a substrate atom of mass M,see Fig. 3, can be written in a Morse 
form 



W(s)=& 



1 _ e a ( s - s b) 



-6, 



with s the internuclear separation, viz. s = F — X , Sb the equilibrium dis- 
tance, 1/a a characteristic interaction length and the binding energy. At 
i = 0, when the emission takes place, the substrate atom receives a sudden 
impulse from other substrate atoms in the collision cascade and the evolution 
of the emitted atom along the normal surface direction may be calculated by 
Energy and Momentum conservation. We have calculated the average velocity 
of the outward motion as function of Vf 



V/ + 



(41) 



in which v 2 = 2 M ^^ ^- Note that in this very simple model other substrate 
atoms do not interact with the QM. For typical binding energies of Alkalis and 
Copper on clean surfaces between 5 and 20 eV we have 10 5 -< v b ^ 10 6 cm/s. 
Since 7" 1 is a characteristic length of the Anderson- Newns interaction, of the 
order of 1 A, e can take values from 10~ 3 to 10 -2 eV at emission energies from 
1 to 10 eV. Consequently 0.01 ^ E\ ;< 0.1 eV. The majority of the metals of 
interest have Fermi wave vectors in the range 0.5 -< kf ^ 2 A _1 and the average 
position of the substrate atom X, in the interaction region, must be of the 
order of a few surface layers, i.e. 1 ^ X ^ 5 A. More quantitatively, in Fig. 5 
we have reported a graphic comparison of Eqs.(38),(40) as functions of e a f 
with the parameters of an alkali sputtering experiment at emission energies of 
5 eV. Despite the two results exhibiting the same dependence on e a f, Eq.(40) 
seems to excessively underestimate the ionized flux for e a f < —0.2 eV, quickly 
approaching the adiabatic regime. The experimentally determined ionization 
coefficient of Cu + sputtering from clean Cu[12] is plotted in Fig.6 as function 
of the emission velocity, Vf. At energies below 10 eV the velocity behavior of 
Eq.(38) is in good agreement with data while at higher energies, above 30 eV, 
the basic result of Eq.( 40) is still valid. 

Inclusion of temperatures, by Eqs.(36) and (37), allows us to conclude that, at 
least in the limit of our model, surface effects are somewhat less sensitive to 
thermal interactions than Anderson-Newns' hopping mechanism. In Fig. 7 we 
have reported the functional dependence of 1Z + on T, for an emission energy 
of the order of 1 eV, a binding energy of 1.2 eV, i.e. e = 2.7 x 10 -3 eV, by 
numerically integrating Eq.(36), the gray line corresponding to the analytical 
solution of Eq.(37). The energy separation e a f has been set to —0.4 eV, typical 
of Cu + /Cu sputtering, and the range of validity of the analytical result (37) 
increases, from T ^ 500 K to T ^ 850 K, with increasing Si (it must be noted 



19 



anyway that when S\ = 0.4 eV or S\ = 0.8 eV the slowness approximation 
can be hardly applied, ry being of the order of |£ a /| _1 ). As shown in Fig.8, 
Eq.(40) varies over several orders of magnitude with increasing temperature 
while Eq.(37) undergoes small variations, within one order of magnitude (see 
also Fig. 7). On the other hand the latter is strongly influenced by changes in 
the characteristic energy e±. We note that, at high temperatures, above 600 K, 
1Z + is mainly determined by the basic Anderson-Newns model, even in the low 
vj-region, while in the interval 200 < T < 500 K, for 0.08 < e 1 < 0.8 eV, QM 
interactions are competitive with the former. It should also be observed that 
the limit of validity of Eq.(39) is not strictly respected below 300 K, which 
would mean that the result form the basic mechanism, reported in the above 
mentioned figure, is over-estimated. On the other hand our result Eq.(36) 
can be applied at all temperatures, provided that resonant charge-exchange 
through the MO is the leading process. Inclusion of thermal interactions attests 
that the two results here discussed are not so dramatically different as the zero 
temperature theory would suggest. In this prospective, it seems an important 
task for future studies to analyze more quantitatively correlations between 
these behaviors, devoting much attention to all the contributions appearing 
in Eq.(26) and, particularly, to the perturbed self-energy of Eq.(30). 

To conclude we observe that finding a power law dependence of 1Z+ on v/, ob- 
served in much experimental work [26,37], is still an open theoretical problem. 
From the scarce contribution found in literature it looks evident that the sole 
resonant interaction of the atomic state with metal electrons is not sufficient 
to obtain such a behavior. In Refs.[4,12] a solution of the form 1Z+ oc v^ has 
been obtained through the introduction of an effective substrate temperature 
into a semiclassical master equation for (n a ), while in Refs. [20,22] the result 
1Z + oc Vf has been proposed as solution of an Anderson-Newns Hamiltonian 
in the screened scattering potential of a moving source localized at atomic 
distances from the surface. These papers have provided the grounds from the 
present analysis, in which non Anderson interaction have been more accu- 
rately modeled by the transient QM. A different approach [25], taking into 
consideration the screening of the nucleus field of a sputtered atom by elec- 
trons simultaneously tunneling in its vicinity through the potential barrier (see 
Figs. 4) , has yield a velocity dependence of the form R + = Hvj (1 + Bvf). In 
the same paper, the author has pointed out that the formula does not hold at 
high emission energies. In the present work, we have explained the different 
behavior of Eqs.(36), (40) and Eqs.(38),(39), in the framework of the same 
model, as depending on the energy received by the sputtered particle in the 
collision cascade. On the contrary, in Ref.[25], the two results are presented 
in relation to the use of two different approaches, the simultaneous resonant 
tunneling and the Anderson-Newns, the second of which is put in doubt be- 
cause of the low v/ features of 1Z+. Our study attests the correctness of this 
latter, in its bare form, providing an implementation that may explain some 
incongruencies with experiments in the critical low Vf region at low- medium 
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temperatures. 



Finally, it is worth to mention that a Hamiltonian resembling Eq.(4), with two 
nearly stationary and directly interacting atomic states, has been used in a few 
different contexts. In Ref . [38] a similar model has been applied to the study of 
resonant charge transfer of scattering proton from an adsorbate-covered metal 
surface. There the authors have modeled different atomic dynamics, regard- 
ing the substrate atom as fixed, thus allowing sensible overlapping between 
the atomic states, and have calculated the ionization coefficient as function of 
the characteristic time of emission, here denoted by ry. The analytical deriva- 
tion has been made on an effective Anderson- Newns Hamiltonian obtained 
by a canonical transformation, thus formally eliminating the second discrete 
state, and compared to simulations. It can be easily shown that the ionization 
probability of Ref. [38], once coupling terms between the discrete states have 
been switched off, reduces to a particular case of Eq.(26), of the presnt pa- 
per, provided that the substrate atom remains fixed, i.e e m ,Vkm (or V m k) and, 
consequently, T™ k , are approximately time- independent. Then by the same 
methods employed in Sec. 4 one can obtain Eq.(37) in the X — > limit. In 
Refs. [39,40] the two state model is used to explain the modification of singlet 
excited states 2 1 S' and 2 1 P of He* interacting with an Al surface. The problem 
is approached through the introduction of an effective 2x2, non Hermitian 
Hamiltonian [41], in the basis of the two, strongly coupled, discrete states, 
in which the non Hermitian terms take into account the interaction with the 
continuum. 

In this paper, we have integrated the reference Hamiltonian of two weakly 
coupled quasi- stationary states directly under the slowness hypothesis and we 
have proposed a manageable expression for the ionization probability that can 
be easily compared with the well known result from the basic theory. Thus we 
have estimated the influence of the new mechanism providing plausible argu- 
ments that substrate excitations produced by the primary ion bombardment 
during sputtering can, in some cases, alter the final charge-state of emitted 
particles. 
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Figure Captions 

• Fig.l: Valence state of a sputtered atom from a metal surface, with time 
dependent energy e a (t) and lorentzian broadening A a (t) , caused by in- 
teraction with an infinitely broad conduction band, with spectrum {et}, 
characterized by a Fermi level £/ and a work function <fi. Surface perturba- 
tions are modeled by the introduction of a quasi molecular orbital, of energy 
e m (t) and broadening A m (t), transiently created during the collision. 

• Fig. 2: Electronic excitations generated by (Fig. 2a) the Anderson-Newns res- 
onant potential, (Fig. 2b) the hoping interaction with the promoted MO and 
(Fig. 2c) scattering from the screened cores of the transient QM. (Fig. 2d) 
first order correlation amplitude between the localized states \a) and \m). 

• Fig. 3: Collision cascade generated by a primary incident atom within the 
first layers of the surface region of a metal substrate. The number of sec- 
ondary emitted particles is mainly determined by binary collisions with sub- 
strate atoms that have received a sudden impulse in the process and form 
Quasi Molecular systems with emitted atoms. The nature of these systems 
is transitory because the recoil velocity of substrate atoms moving within 
the metal is larger than the emission velocity. 

• Fig. 4: Schematic representation of the effective one-electron potential (black 
line) of the first quantization Hamiltonian (2) as a function of the normal 
position coordinate to the surface plane. Three main contributions can be 
distinguished taking into account the effect of surface barrier v s , the sput- 
tered atom v a , and the quasi-molecule v m . Note that as the interatomic 
separation in the QM increases, the state \m) approaches the continuum 
losing localization and the interacting terms Vk m and Vkw become negligi- 
ble. 

• Fig. 5: Comparison of Eqs.(38), black line, and (40), gray line, as functions of 
the energy separation e a f = e a — Ef for a typical alkali sputtering experiment 
at E = 5 eV of emission energy and v b = 7.810 6 cm/s of binding velocity. 
The average outward velocity has been evaluated from Eq.(41), the interac- 
tion length 7 _1 set to 0.5 A , the product kfX to 5 and the characteristic 
energy e\ to 10eo- 

• Fig. 6: Experimental data of Ref. [12] fitted, in the low velocity range vj -< 
6 x 10 6 cm/s, with a linear dependence of 1Z + on v (black line). The binding 
velocity vj, has been used as a fitting parameter. Its result, in agreement 
with expected values, has been inserted into Eq.(39) and used to fit data in 
the range Vf y 9 x 10 6 cm/s (gray line). 

• Fig. 7: 7£ + -dependence on temperatures in the case of dominant charge- 
exchange processes through the MO, for different values of the inverse of 
the average life-time of the QM, e±. In each figure the black line has been 
obtained by numerically integrating Eq.(36) and inserting the result into 
Eq.(31) and the gray line represent the low temperature behavior of Eq.(37). 
eo has been fixed for an emission energy of 1 eV, with the v& of Fig. 5 and e a f 
to the literally value[12] —0.4 eV. Despite a weak variation with increasing 
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temperature, 1Z+ is strongly sensitive to S\ and the range of validity of the 
analytical solution (37) increases with increasing e±. 
• Fig. 8: Comparison of different temperature behaviors of the ionization prob- 
ability in the two limiting cases discusses in the article: The gray line corre- 
sponds to the numerical solution for 1Z + , in absence of surface interactions, 
see Eq.(40), and the black lines have been obtained as in Fig.7. Parameters 
of Fig.7 have been used. Note that the Anderson- Newns hopping mechanism 
seems to prevail at high temperatures when T y 600 K. 
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